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\l^ , We consider one scenario of transformation of a cycle of a differential equa- 

te ' tion, which we call the blow up. In this scenario, a cycle born via the Hopf 

bifurcation grows to infinity as a parameter ranges over a finite segment. 
From another perspective, in the product of the phase space and the pa- 
rameter axis there is a branch of cycles connecting the equilibrium and 
infinity. We first discuss the existence of such a branch for planar differen- 
1^^ ' tial equations where the proof can be based on the Poincare theorem. A 

jrt ' Lotka-Volterra type system with a competition-cooperation term is consid- 

ered as an example. Then we discuss the existence of a branch of cycles 
stretching from zero to infinity for a class of higher order quasilinear equa- 
tions: this theorem continues the results of Refs. 1,2. The results for planar 
systems are presented in the next section. Section 3 contains the main result 
for higher order equations. We briefiy sketch some points of the proofs. 



o 
o 



1. Introduction 



2. Planar systems 

Consider a planar system 

x' = f(x,y;X), y'^g{x,y;X) (1) 
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with a parameter A G [0, 1], where the functions f,g are continuously dif- 
ferentiable with respect to the phase variables x,y and continuous with 
respect to the set of all arguments. We say that the system has a branch of 
cycles connecting the origin and infinity if for any bounded open domain 
G 9 (0,0) of the phase plane the system has a cycle (for some A G [0, 1]) 
that belongs to the closure G of G and touches the boundary dG of G. This 
definition, close to the classical weak definition of continuous branches of 
fixed points'^ , is discussed in Rcf. 1. 

Proposition 2.1. Suppose that the following conditions hold: 

(i) /(0,0; A) == and the origin is the only equilibrium point of system (1) 

for all \ e [0,1],- 
(a) the Jacobi matrix J = J(A) of system (1) at the origin is invertible for 

allXe [0,1]; 
(Hi) the spectrum of J is in the open right half-plane for A = and in the 

open left half-plane for A = 1; 
(iv) the system does not have cycles for A = and A = 1. 

Then system (1) has a branch of cycles connecting the origin and infinity. 

Condition (ii) ensures that no equilibrium branches from the origin, 
which agrees with condition (i). According to condition (iii), the origin 
changes stability as A ranges over the segment [0, 1]. Conditions (ii) and (iii) 
guarantee that the Andronov - Hopf bifurcation occurs within this segment. 
Condition (iv) is satisfied, for example, if the system has a global Lyapunov 
function for A = and A = 1. Under the assumptions of Proposition 2.1 
the system does not have homoclinic orbits; the cycles surround the origin. 

As an example, consider the Lotka-Voltcrra type equations 

x' ^x{a-by), y' = y{-c + dx + f{y; X)); x,y > (2) 

with positive parameters a, b, c, d. The last term / in the predator equation 
accounts for competition (/ < 0) or cooperation (/ > 0) in the predator 
population. We assume that the behavior of the predator can depend on 
its number: it competes if the population number is above a certain thresh- 
old y — y{X) and starts to cooperate when the population falls below the 
threshold, hence /(y; A) < for y > j/(A) and f{y; A) > for y < y{X). For 
example, / = arctany — Ay with A € (0, 1). 

Suppose that c > f{ab^^; A) for all A (the condition c > f means that 
the predator extincts in the absence of prey). Then system (2) has a unique 
positive equilibrium (x*,y,) = {d^^{c— f{ab~^;X)),ab~^) and Proposition 



July 27, 2010 1:0 WSPC - Proceedings Trim Size: 9in x 6in rachinskii'dmitrii 



2.1 can be applied after the logarithmic coordinate transformation and the 
shift of the equilibrium to zero. The equilibrium (x*, j/*) is stable if /^ < 
and unstable if /^ > at y = y» , where fy — df/dy is the partial derivative 
of ./ = f{y. A). Hence conditions (i) - (iii) of Proposition 2.1 are satisfied if 

/;(y,;0)>0, /;(y,;l)<0. (3) 

If in addition 

(y - y*)ifiy; o) - ./(y*; o)) >o,y>o, y^y,, (4) 

(2/ - y*)(/(2/; 1) - fiy*: l))<0,y>0, y^ y,, (5) 

then condition (iv) is also satisfied, because in this case V ~ (x — x^, liix)d+ 
{y — y* lny)6 is a Lyapunov function of system (2) for A = 0, 1 with V ~ 
[y — y*){f{y; A) — /(j/*; A))6. Hence, relations (3) - (5) ensure the existence 
of a branch of cycles connecting the equilibrium and infinity for system 
(2) (the cycles lie in the positive quadrant x,y > where the system is 
defined). In particular, these relations hold if / strictly increases for A = 
and strictly decreases for A = 1, as in the above example / = arctany — Ay. 
Numerical simulations confirm that a stable positive cycle born via the 
Hopf bifurcation blows up to infinity for this / and demonstrate the same 
effect for other competition-cooperation terms, such as / = y — Xy^ or 
f — y^ — Ay"', included in the equation for predator, or prey, or both. 

To prove Proposition 2.1, one can first note that under its conditions 
the equilibrium can not have eigenvalues of different sign and consequently 
the system does not have homoclinic orbits. Using the Poincare theorem, 
one derives from this fact that if the system has an orbit 7 in a bounded 
domain G with 7 n dG ^ 0, then it also has a cycle C (iG with C fl dG 7^ 0. 
To complete the proof by contradiction, assume that there is no such a 
cycle and hence no such an orbit 7 for some bounded domain G 9 0. 
Consequently, S\ fl dG = 0, where S\ denotes the invariant set of the 
system in the domain G (note that G S\). Therefore G is an isolating 
neighborhood for S\ and the Conley index Ind S\ of S\ with respect to G 
is defined^ . Because the system has no homoclinic orbits and, by condition 
(iv), there is no cycles for A = 0, 1, the Poincare theorem implies that 
Sq = S\ ~ 0. Moreover, condition (iii) implies Ind 6*0 7^ Ind 6*1. This, 
however, contradicts the invariance of the Conley index under homotopic 
transformation of the vector field: Ind S\ should be the same for all A for 
any isolating neighbprhood G oi S\. Given any open bounded G 9 0, this 
contradiction proves the existence of a cycle C C G with C n dG ^ for 
some A, i.e. the conclusion of the proposition. 
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3. Quasi-linear higher order equations 

Consider the equation 

i(|;A)x-/(x;A), (6) 

where L{p;X) = p^ + ai{X)p^^^ + • • • + ao(A) is a polynomial with con- 
tinuously differentiable coefficients. Assume that the continuous function 
/(x; A) satisfies /(O; A) = and the global Lipschitz estimates 

\f{xi;X)~f{x2;X)\<k\x,~X2l \fix;Xi)-fix;X2)\<l\x\\Xi-X2\; (7) 

hence the equation has the zero solution a; = for all A. Define the matrix 

, , / ^L'^{X;wi) — ^m Lp{X; wi) 
^'^' '^\^^L'^{X;wi) ^L'piX;wi) 

We say that equation (6) has a Lipschitz continuous branch of cycles con- 
necting zero and infinity if there are Lipschitz continuous functions X(r), 
w{r) with values in segments [A_, A+J, [w-, w+] (w^ > 0) such that for ev- 
ery r > equation (6) with A = A(r) has a periodic solution Xr{t) ~ x{t; r) 
of the period 2T:/w{r), the function x{t; r) is Lipschitz continuous in r and 

||2;r||c~^0 as r — > 0, ||xr||c ~^ oo as r — ;■ oo. 

Theorem 3.1. Assume that for some q > the relation \L{wi; X)\ < q 

defines a simply connected bounded domain Dq on the plane {w,X), the 
equation L(wi; X) = has a unique solution {wq,Xo) in Dq, the matrix 
J(w, A) is nondegenerate in Dq, and L{nwi; X) =/= in Dq for any integer 
n 7^ ±1. Then there are sufficiently small k,l > Q such that equation (6) 
with any function f satisfying the estimates (7) has a Lipschitz continuous 
branch of cycles connecting zero and infinity. 

The method of the proof of Theorem 3.1 leads to explicit estimates of 
the Lipschitz coefficients k,l, which ensure the existence of the branch of 
cycles connecting zero and infinity. 

A natural parameter r is the amplitude of the first harmonics of the 
periodic solution. Theorem 3.1 can be proved by contraction mapping prin- 
ciple. To construct the corresponding mapping, let us first note that for 
any (w,A) G Dq the differential operator L{w-^;X) with the 27r-periodic 
boundary conditions is invcrtiblc on the codimension 2 subspace E of 
l7 = L^(0, 27r) which is orthogonal to sint and cosi (this operator, how- 
ever, is not invertible on the whole space L^, because L{woi',Xo) — 0). 
Secondly, the planar map (w. A) i— > {^L{wi; X),'<smL{wi; X)) is invertible 
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on Dq. Now, denote by Pg and Pc the orthogonal projectors onto sint and 
cos i in L^ , define the orthogonal projector Q = I — Pg — Pc onto E = QL^ 
and consider the space of triples {u,v,y) G M x M x QL^ with the norm 



||(u,w,2/)|| = \/ u^ + v'^ + \\y\\L2 , where y ~ j/(i). In this space, for each 
value of the parameter r > consider the mapping 

(u, V, y) ^ Ar{u, V, y) = r-\Pg,f{x(t); X),PJ{x{t); A), Q/(x(t); A)), 

where x — x(t) and A are defined by the relations 

M = 3ffeL(wi;A), v ^ ^rn L{wi; X) , x(i) = r(7r"^/^ sini + /i(t)), (8) 

L{w^]\\h{t) = y{t), h{0)^h{2T:),h'{0)^h'{2T:), heQh^. (9) 

Due to the invertibility of these relations mentioned above, the mapping 
Ar is well-defined for all u^ + v'^ < q^ and all y G QL^. The definition of Ar 
ensures that every fixed point of Ar defines a 27r-periodic solution x = x{t) 
of the equation L{w-^;X)x — /(x; A) by the formulas (8), (9) and hence a 
27r/'u;-periodic solution x{wt) of equation (6). The proof is completed by 
showing that if fc, I are sufficiently small then Ar is a contraction on the 
ball \\{u, V, y)\\ < q, which is invariant for Ar, for each r > 0. 
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